Long-range spin-triplet proximity effect in Josephson junctions with multilayered 

ferromagnets 
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We study the proximity effect in SF'(AF)F'S and SF'(F)F'S planar junctions, where S is a clean 
conventional (s-wave) superconductor, while F' and middle layers are clean or moderately diffusive 
ferromagnets. Middle layers consist of two equal ferromagnets with antiparallel (AF) or parallel (F) 
magnetizations that are not coUinear with magnetizations in the neighboring F' layers. We use fully 
self-consistent numerical solutions of the Eilenberger equations to calculate the superconducting pair 
amplitudes and the Josephson current for arbitrary thickness of ferromagnetic layers and the angle 
between in-plane magnetisations. For moderate disorder in ferromagnets the triplet proximity effect 
is practically the same for AF and F structures, like in the dirty limit. Triplet Josephson current is 
dominant for d! ~ hvp /2h' , where d' is the F' layer thickness and h' is the exchange energy. Our 
results are in a qualitative agreement with the recent experimental observations [T. S. Khaire, M. 
A. Khasawneh, W. P. Pratt, and N. O. Birge, Phys. Rev. Lett. 104, 137002 (2010)]. 

PACS numbers: 74.45.-l-c, 74.50.-|-r 



In hybrid systems containing superconducting and fer- 
romagnetic metals, triplet correlations are inducedji"— In 
the case of homogeneous magnetization, the triplet pair 
amplitude has zero total spin projection on the magne- 
tization axis: Fto{t - t') oc|t (*) I (i'))+ 14- (^') t (0)- 
For equal time t = t' , this function vanishes, in agree- 
ment with the Pauli principle. Therefore, Fto is an odd 
function of the time difference or equivalently, odd in fre- 
quency. The exchange field mixes the usual spin singlet 
pairing correlations, Fg oc|t4,)— |4,t), and the spin triplet 
Fto because the wave functions |t4-) and |4,t) acquire rel- 
ative phase difference in the ferromagnet.- 

Amplitudes Fg and Fto are short-ranged and oscillate 
spatially in the ferromagnet, both in the clean and the 
dirty limit. '^ In the clean limit, they decay algebraically 
in the ferromagnet with hvp/h, where h is the exchange 
energy and vp is the Fermi velocity. In the dirty limit, 
they decay exponentially with the characteristic length 
y/hD/h, where D = vpl/3 is the diffusion coefficient, and 
I is the electron mean free path. However, this is not true 
for the clean single-channel junctions where all pairing 
correlations are long-ranged and their spatial decay in 
the ferromagnet is independent of the exchange fieldfi 

The situation is quite different for inhomogeneous 
magnetization: not only Fs and Fto amplitudes exist 
but also odd triplet pair amplitude Fti with ±1 to- 
tal spin projection on the magnetization axes emerges 
in the ferromagnet ii"— The triplet component Fti is 
not suppressed by exchange interaction and penetrates 
into the ferromagnet over large distance on the order of 
y^hD/ksT in the dirty limit (/ < hvp/h) and moderately 
disordered ferromagnets, and, likewise, over the distance 
hvp/kBT in the clean limit {I > hvp/ksT). 

It is not difficult to understand why Fti triplet com- 
ponent is induced in inhomogeneous ferromagnets. Con- 
sider a system where inhomogeneity is represented by 
two ferromagnetic layers with noncoUinear magnetiza- 
tions and angle a between them. As we have already 
pointed out, in each ferromagnetic layer the exchange 



field generates Fto from Fg. When Fto correlation pene- 
trates into neighboring ferromagnetic layer it mixes with 
Fti having non-zero total spin projection due to differ- 
ent orientation of magnetizations. This implies that Fti 
triplet component is generated in the presence of inhomo- 
geneous magnetization and is proportional to Fto sin a at 
the interface between two ferromagnetic layers. There- 
fore, for fully developed triplet proximity effect one of 
two layers should be sufficiently thin to provide large Fto 
at the interface between the layersi^ 

Another possibility for dominant long-range triplet 
component was suggested in Refs. ^-_10. In that ap- 
proach the role of thin ferromagnetic layers is replaced 
by spin-active interfaces described by scattering matrix. 
The elements of the scattering matrix can be considered 
as phenomenological parameters. Purely microscopic 
approach that retains quasiparticle information at the 
atomic scale, with spin-dependent scattering potentials 
at the interfaces, was considered in Ref. [ill 

Besides the impact on the Josephson current, another 
signature of odd-frequency triplet pair correlations has 
been proposed recently: The density of states in the 
presence of the odd-frequency pairing is enhanced, and 
acquires a zero-energy peakj ^^d^ 

Experimental results that may support the existence 
of long range triplet amplitudes have already been 
obtained. ^"^"^^ However, in these experiments it was not 
possible to tune the junction parameters, and the tran- 
sition from usual singlet to long-range triplet proximity 
effect has not been observed. 

Quite recently, long-range triplet proximity effect has 
been observed in SF'(AF)F'S junctions with F' layer 
thickness as a controllable parameter?!! Here F' is a 
weak ferromagnetic layer (PdNi or CuNi) and AF is 
synthetic antiferromagnet: an exchange-coupled trilayer 
Co/Ru/Co. Previously, a similar arrangement (with ho- 
mogeneous middle layer F) has been analyzed theoreti- 
cally and proposed as a good candidate for the long-range 
triplet effect ,^ii^ 



Junctions with F and AF middle layers were analyzed 
recently for the case when F' layer thickness is much 
smaller than ^JhD/h' }- However, within this approxi- 
mation the triplet component is not dominant and con- 
sequently results are not applicable to the experiment 
(Ref. 17). More interesting case was considered by the 
same authors, when F' layer thickness is arbitrary but 
the middle layer is homogeneous, as in Ref. J^. In the 
dirty limit (and moderately disordered ferromagnets as 
we will show), results are practically the same for F and 
AF structures of the middle layer. The situation is dif- 
ferent for the case of clean ferromagnets. 

In this paper, we study the proximity effect in clean 
or moderately diffusive SF'(AF)F'S and SF'(F)F'S pla- 
nar junctions, where S is a conventional (s-wave) super- 
conductor, while F' and middle layers are ferromagnets. 
Middle layer consists of two equal ferromagnets with an- 
tiparallel (AF) or parallel (F) magnetizations that are 
not coUinear with magnetizations in the neighboring F' 
layers. 

To calculate the Josephson current and pair correla- 
tions in the clean limit and for moderately diffusive fer- 
romagnets, we use the Eilenberger equations^ for a junc- 
tion along the a;- axis 



hvxdxg + uinh - iV + K + h{g)/2T, g 



0, 



(1) 



with normalization condition g^ ^ \. Disorder is charac- 
terized by the average time r = l/vp between scattering 
on impurities. We indicate by • " • and • " • 2x2 and 4x4 
matrices, respectively. Here, Vx = vp cos 9 where 9 is an 
angle between the Fermi velocity and the x-axis, fi are 
the Pauli matrices in the particle-hole space, the brack- 
ets (. . .) denote angular averaging over the Fermi surface 
(integration over 0), and [ , ] denotes a commutator. The 
quasiclassical Green functions are given by 



gs +gt • o- ifs + ft- cr)icTy 



(2) 



where <t = {a'x,d'y,az) are the Pauli matrices in spin 
space. With this parametrization,^ it is clear that gs and 
fs remain unchanged under rotation of magnetizations, 
while Qt and ft transform like ordinary vectors. The 
conjugate Green's functions satisfy /s(i^n) = fl{~^n) 
and ftioon) = —ft{~^n)-'^^ The exchange field term in 
ferromagnets is given by F = Y{,e[h{x)-a-]+iT^ lTa[h{x)-a-] 
and the pair potential in superconductors is A = (t+ A -f- 
f_ A*)(7jy, where t± = Tx ^ ify. The exchange field h has 
the following x-dependence 



h 



h'{— sin a y + cos a z), 

h z, 

±h z, 

h'{sm a y + cos a z). 



—d ~ d! < X < —d, 
-d <x <0, 
< X < d, 
d< X < d + d', 



(3) 

where d' and 2d are F' and (F) or (AF) thickness, re- 
spectively. The angle between magnetizations in F' and 
neighboring ferromagnetic layers is a; the plus (minus) 
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FIG. 1. (color online) Spatial dependence of singlet and 
triplet pair amplitudes Fs, F^ and F^, normalized to the 
bulk singlet amplitude F^b, for T /Tc = 0.1, h' /Ef = 0.05, 
h/Ep = 0.1, 2d = SOOifc^S d! = 20fc^\ a = 7r/2 and two 
values of the mean free path: I = oo (dashed curves) and 
I = 200fc]^^ (solid curves). All amplitudes are calculated for 
the ground state, (j) = Q. The SF'(AF)F'S junction geometry 
is shown in the background; arrows and circles show orienta- 
tion of magnetizations in each layer. 



sign is for F (AF) middle layer. In the absence of out-of- 
plane magnetization, the amplitude {ft)x — 0. 

The supercurrent flowing through the junction is given 
by the normal Green function 



/((/.) = 'KeN{Q)SkBTlniY^{vxgs{vx)). 



(4) 



where <j) is the macroscopic phase difference across the 
junction, iV(0) is the density of states per spin at the 
Fermi surface, and S is the area of the junction. Equal- 
time pair amplitudes are defined in terms of anomalous 
Green functions as 



F, = 



iv(o)fcBr^(/,), 

F< = T,N{Q)kBT J2 ((/t)^>> 

w„<0 

NiO)kBT J2 iift)y)- 



Ft-i 



(5) 
(6) 
(7) 



w„<0 



Equal-time triplet amplitudes identically vanish accord- 
ing to the Pauli principle, hence we defined auxiliary 
functions using summation over negative frequencies 
only. The spatial variation of the time-dependent triplet 
pair amplitudes is qualitatively the same as for auxiliary 
functions. 



We consider the case of fully transparent interfaces, 
i.e., strong proximity effect. The opposite case of low 
transparency was considered in Ref. [la- Using continuity 
of the Green functions at interfaces, Eqs. ([T]) are solved 
iteratively with the self-consistency condition A = XFg. 
Iterations are performed until self-consistency is reached, 
starting from the stepwise approximation for the pair po- 
tential A = A(r)[e-*^/2e(-x -d-d') + e''l>/^Q{x -d- 
d')]. The temperature dependence of the bulk pair po- 
tential is given by A(T) = A(0)tanh (l.74:y/Tc/T - iV 

For arbitrary mean free path in ferromagnets we employ 
the iterative procedure starting from the clean limit. 

Figure [T] shows the spatial variation of the pair ampli- 
tudes for the SF'(AF)F'S junction geometry with magne- 
tizations in F' layers orthogonal to the neighboring mid- 
dle layers, a — 7t/2. In this case the 0-state is the ground 
state. Transition to the 7r-state can be tuned with rela- 
tive orientation of magnetizations in F' layers (7r-state is 
the ground state for parallel magnetizations in the two 
F' layers). For corresponding SF'(F)F'S junctions pair 
amplitudes are practically the same. For (j) — 0, the sin- 
glet and the long-range triplet amplitudes, Fs and Fti, 
are real while the short-range triplet amplitude Fto is 
imaginary. 

The F' layer thickness d' = 20k p^ is chosen to give the 
maximum triplet current for moderate disorder in ferro- 
magnets {I = 200k p^). This explains the large difference 
between amplitudes of the long-range triplet component 
in the clean and moderately disordered case. Fig. [T] 

We observe oscillatory decay of Fg and F^ ampli- 
tudes, dependent on the exchange field. In contrast, long 
range triplet component, F^, is monotonic in the mid- 
dle layer and its decay length is independent of the ex- 
change field magnitude. Thin F' layers are considered 
as weak ferromagnets, h' / Ep = 0.05, and thick middle 
layers, 2d = 500fc^^, represent strong ferromagnets. For 
illustration of pair amplitudes, h/Ep =0.1 is used in or- 
der to have larger period of spatial oscillations, although 
h/Ep = 0.3 is used in other illustrations. All ampli- 
tudes are normalized to the bulk singlet pair amplitude 
F,, = 2nN{0)kBTj:^^ A/^u^ + A^. 

Next we examine the Josephson critical current de- 
pendence on F' layer thickness. This quantity is actu- 
ally measured in the experiment J^ Figure [5] illustrates 
Ic{d'), normalized to the resistance i?jv = 2TT'^h/Se^kp, 
for three types of junctions: SF'(AF)F'S, SF'(F)F'S and 
SNS, where N is the corresponding normal nonmagnetic 
metal {h — h' — 0). Here, we take mean free path 
I = 200kJ, in ferromagnetic and N metals. 

In the clean limit, Ic{d') curves for SF'(AF)F'S and 
the corresponding SNS junctions coincide,^ while in the 
dirty limit Ic{d') curves for SF'(AF)F'S and SF'(F)F'S 
curves are practically the same.^^ For intermediate disor- 
der in ferromagnets, the critical current is always larger 
in SF'(AF)F'S than SF'(F)F'S junctions. Fig. H 

We find for moderate disorder in ferromagnets that 
the largest critical current is almost as big as for the cor- 
responding SNS junction. Note that in the dirty limit 
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FIG. 2. (color online) Dependence of the Josephson critical 
current on the thickness d! of F' ferromagnetic layer , for 



500fc^ 



T/Tc = 0.1, h'/Ep 



2d 

200fc^^, and for three types of junctions 

and SF'(F)F'S. 
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FIG. 3. Dependence of the Josephson critical current on 



the middle layer thickness 2d, for d' 



20fc; 



T/T, = 0.1, 



h'/Ep = 0.05, h/Ep = 0.3, I = 200fc^\ and for two values of 
the misorientation angle: a = 7r/2 (thick and thin solid lines 
for AF and F structures) and for a = (AF structure, dashed 
line). Arrows and circles show orientation of magnetizations 
in ferromagnetic layers. 



maximum triplet critical current is only 12% of the crit- 
ical current of the corresponding SNS junction.^* The 
position of maxima of Ic{d') strongly depends on /i', 
'^'max ~ hvp /2h' . It is practically independent of h and 
d, in agreement with the results of Ref. |l^. 

The critical Josephson current dependence on the mid- 
dle layer thickness is shown in Fig. |3l The F' layer thick- 
ness is set to d'^aax = 20k p^. We consider Ic dependence 
on 2d for two values of misorientation angle: a — 7r/2 
(fully developed triplet proximity effect, solid lines) and 
a = (no long range triplet component, dashed line). 
For thick middle layer, 2d = 600fc^\ in both AF and F 
geometries, Ic is an order of magnitude larger when the 
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FIG. 4. The current-phase relation /(<^) for SF'(AF)F'S junc- 
tion: d! = 20kp\ 2d = 600kp\ T/n = 0.1, h'/Ep = 0.05, 
h/Ep = 0.3, I = 200kp^, and a = n/2. 



triplet component Fti is present. These results, Figs. [2] 
and[3l are in a qualitative agreement with experimental 
observatioufii 

The current-phase relation is almost harmonic for 
a moderately disordered SF'(AF)F'S junction {I = 
200kp^), Fig. m Here, the F' layer thickness is optimal 
for long-range triplet Josephson effect; i.e., the usual sin- 
glet Josephson critical current is an order of magnitude 
smaller. For the in-plane magnetizations considered here 
we did not find any unusual I{(f)) dependence, in agree- 
ment with Refs. ,6, and ,18,, 

However, we expect unusual /(0) dependence for the 
samples used in the experiment (Ref. llTl ). Supercurrent 



count the possibility of inhomogeneous out-of-plane mag- 
netisations in PdNi layers and the fact that transport 
properties of the minority and majority electrons at the 
Fermi surface in Co layers are very different. Quasi- 
classical approximation we used in this paper does not 
allow for the latter possibility; this is the principal limi- 
tation of our approach. 

A new long-range triplet component, {ft)x — (/tt ~ 
/4.4.)/2, is generated in the case of non-zero relative lon- 
gitude angle x between magnetizations in two F' layers 
(i.e., inhomogeneous out-of-plane magnetization). This 
can be readily seen from Eq. ([T]) for {St)x component. 
The presence of {^ft)x implies different triplet Josephson 
current flow for majority (J-^) and minority (/j,) electrons. 
In the circuit theory approximatiouji^ it was obtained 
that /^_^ = C^^^sin ((^ ± x), where C^^^ are proportional 
to densities of states and Fermi velocities of spin sub- 
bands. Although the case of very different Fermi veloc- 
ities of the subbands is not accessible within the circuit 
theory approximation, it is reasonable to assume that 
similar current-phase relations are valid. Consequently, 
there is a finite supercurrent at zero phase difference, as 
was predicted for the half-metallic middle layer.— 

It would be very interesting to measure the /((/)) de- 
pendence for the samples used in the experiment (Ref. 
[l7r ). since the out-of-plane magnetizations in thin PdNi 
layers are their typical feature. ^^ A non-zero supercur- 
rent for (/) = could provide an unambiguous evidence of 
the triplet proximity effect. 

We acknowledge useful discussions with Norman Birge, 
Marco Aprili, Ivana Petkovic, Mihajlo Vanevic and Ivan 
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istry of Science, Project No. 141014. 
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